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We present preliminary results for the AI = 3/2 kaon decay matrix elements using domain wall 
fermions and the DBW2 gauge action at one coarse lattice spacing corresponding to a -1 = 1.3 
GeV. We employ an extention of the Lellouch and Luscher formula for non-zero total momentum 
to extract the infinite volume, center-of-mass system decay amplitudes. We compare the result 
of ReA2 with previous results calculated by several methods, and experiment. We also show the 
1 = 2 %n scattering phase shift. 
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1. Introduction 

Non-leptonic kaon decay includes interesting subjects such as, the AI = 1/2 selection rule and 
the CP violation parameter e'/e. However, it is difficult to calculate the K —> Tin weak matrix ele- 
ment directly on lattice due to difficulty of calculation of the two-pion state in finite volume, which 
was pointed out by Maiani and Testa [p. So far many groups employ an indirect method [||] using 
Chiral perturbation theory (ChPT) to avoid the difficulty In the indirect method K — > %% process 
is reduced to K — > % and K — ► processes. Since in the decay process the final state interaction of 
the two-pion is expected to play an important role, the reduction may cause systematic errors. 

For the direct method, where the two-pion state is calculated on lattice, we need complicated 
calculations and analyzes, e.g., diagonalization of a matrix of correlation functions [gj], to treat the 
two-pion state with non-zero relative momentum on lattice. However, recently Kim [Q] reported an 
exploratory study of the direct calculation with H-parity boundary conditions, where complicated 
analyzes are not required. He succeeded to extract the two-pion state with non-zero relative mo- 
mentum from the ground state contribution of correlation functions, because the two-pion state with 
zero momentum is prohibited by the boundary condition. In the non-zero total momentum (Lab) 
system, |P| / 0, the ground state of the two-pions is \%(0)ii{P)), which is related to the two-pion 
state with the non-zero relative momentum in center-of-mass (CM) system. Thus, we can extract 
the two-pion state with non-zero momentum from the ground state contributions as well as in the 
H-parity boundary case. However, we cannot naively apply the formula proposed by Lellouch and 
Liischer (LL) ^ to the calculation of the Lab system, because the original formula is derived in the 
CM system. Before the simulation, we have to modify the LL formula to connect the finite volume 
decay amplitude in the Lab system to that of the CM system in infinite volume. Very recently, Kim 
et al. [||] and Christ et al. |Q] suggested the formula which is an extension of the LL formula for 
the Lab system calculation. Here we attempt to apply this extended formula to the calculation of 
the AI = 3/2 kaon weak matrix elements with domain wall fermions and the DBW2 gauge action 
at a single coarse lattice spacing. 



2. Methods 

Two groups [||, [7|] proposed a formula between a decay amplitude in infinite \A | (CM system) 
and finite volume \M\ (Lab system), 

where 7 is a boost factor, E n% is the two-pion energy in the CM system, p 2 = E 2 n /A — m 2 , and 
8 is the scattering phase shift of the final state interaction. The function (j>p with P being the total 
momentum, derived by Rummukainen and Gottlieb is defined by 

tanfe(g) = - Jlf'l (2.2) 

where q 2 = (pL/ln) 2 , and 

4(i;<? 2 ;r) = -±= V 2 , 2 , * , - — ? , (2.3) 
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in the P = (0,0,27l/L) case. The formula eq. (2.1) is valid only for on-shell decay amplitude, i.e. 
E n% = m-K as is LL formula [g]. 

We calculate the four-point function for AI = 3/2 K — > %% decay with total momentum P = 
and (0,0,27r/L). The four-point function G,(?) is denned by 

Gi{t) = {Q\n + Ti-(P,t 7l )o] l2 (t)K°(P,t K )\0), (2.4) 

3 /2 

where the operators O] are lattice operators entering AI = 3/2 weak decays 

0^88 = $°d a ) L [(u b u b ) L , R - [d b d b ) L . R ] + (?'u«) L (u b d b ) Lfi (2.5) 

0l£s = ^'d h )L [{u b u a ) R - (rfV) R ] + (ru b ) L (u b d a ) R (2.6) 

3/2 3/2 

with (qq)L,R = <7%i(l ± Ts)*? and a, b being color indices. 2 j and O gg are the operators in the 

3 /2 3 /2 

(27,1) and (8,8) representations of SU (3)l (g) S£/ (3)r with 1 = 3/2, respectively. O m ' 8g equals 8 g 
with its color summation changed to cross the two currents. We also calculate four-point function 
for two-pions and the two-point function of a kaon with zero and non-zero total momenta, to obtain 
those energies and amplitudes. 

We employ the domain wall fermion action with the domain wall height M = 1.8, the fifth 
dimension length L s = 12 and the DBW2 gauge action with j3 = 0.87 corresponding to a -1 = 1.3 
GeV. The lattice size is L 3 • T = 16 3 • 32, where the spatial volume corresponds to about 2.4 fm. 
We use four light quark masses, m u = 0.015,0.03,0.04 and 0.05, for the chiral extrapolation of 
the decay amplitudes, and six strange quark masses, m s = 0.12,0.18,0.24,0.28,0.35 and 0.44, for 
the interpolation of the amplitudes to the on-shell point. We fix the two-pion operator at t % = 0, 
while we employ three source points = 16,20 and 25 for the kaon operator to investigate the 
dependence of the statistical error of the Lab system decay amplitude and to check the consistency 
of these results. The numbers of configuration used are 1 1 1 for tjc = 16 and 20, and 100 for t% = 25. 
A quark propagator is calculated by averaging quark propagators with periodic and anti-periodic 
boundary conditions for the time direction to obtain a propagator with 2T periodicity. We also use 
a wall source with Coulomb gauge fixing for the quark propagators. 

3. Results 

In order to utilize the eq. ( |Q1 ) we need to measure derivatives of the scattering phase shift for 



the final state and evaluate the function (f)p(q) denned in eq. (2.3). The phase shift 8(p) is evaluated 



from the two-pion energy by the finite volume formula We carry out a global fitting of 

tan5(p) E„„ 

T{m 7t ,p) = — (3.1) 

P 2 

with a polynomial function a\Qin 2 n + fl , 20 w 4 + a 0\P 2 t0 obtain the derivative of 8 with respect to p. 
The measured value of T (left graph) and 8 (right graph), and their fitting results are plotted in 
Fig. [T[ The prediction of 8 from ChPT with experiment is also shown in the right graph. 

To determine off-shell decay amplitudes in finite volume, we define ratio of correlation func- 
tions Ri(t) for i = 27, 88 and m88 as 

V3Gi(t)Z KJ[ Z K 

Riit) ~ GUt)G K {t) ' (3 " 2) 



351/3 



AI = 3/2 kaon weak matrix elements with non-zero total momentum 



Takeshi Yamazaki 



T(m/,p 2 )=tan8(p)EJ2p 
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Figure 1: Measured values of 7\ defined in eq. (3.1 ), and 8 obtained from CM and Lab calculations. Solid 
lines are fitting results. 
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Figure 2: Rzj(t), defined in eq. (3.2), with different tx obtained from CM and Lab system calculations. 
Small figures show averaged values in flat regions. 



where G %n {t) and Gx{t) are I = 2 two-pion four-point function and kaon two-point function, re- 
spectively, and Z %% and Zx are the overlap of the relevant operator with each state. The factor yfb 
comes from changing basis to the isospin basis. When these correlation functions are dominated 
by each ground state, the ratio /?/(/) will be a constant for those values of t. We then determine the 
off-shell amplitude in this region. 

The figure ^| shows Rn{t) for the lightest pion mass and fixed strange quark mass obtained 
from the CM and Lab system calculations. We plot the ratios with three different kaon source 
points tx in the figure. In the CM system, the ratios with different tx are reasonably flat and 
consistent with each other in a region t n <C t <C tx- The off-shell amplitudes are determined from 
averaged values in the flat region, whose values are presented in small figure. The three averaged 
values are also consistent with each other. In the Lab system case, the ratios are noisier than those 
in the CM case. While the error with tx = 25 is very large, the error decreases as t% decreases. 
Thus, we choose the results with tx = 16 in the following analysis, because the result is consistent 
with the others, and has the smallest error of the three cases. We determine the on-shell amplitude 
in finite volume by a linear interpolation of the off-shell amplitude with different strange quark 
masses at fixed light quark mass. The decay amplitudes in infinite volume are obtained for the 
extended formula of eq. ( |2.1[ ) by combining the on-shell amplitude and the derivatives of 8 and (j)p. 

The matching factors to determine the weak matrix elements were previously calculated in 
Ref. [Q] with the renormalization independent (RI) scheme by a non-perturbative method at the 
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Figure 3: Weak matrix elements obtained from CM and Lab system calculations. Squre symbols are fitting 
results with p = 0. 
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Figure 4: Measured ReA2 and result from previous works. Square symbols denote fitting results. 



scale /i = 1.44 GeV. The weak matrix elements \Af l \ for i = 27, 88 and m88 in the RI scheme are 
plotted in Fig. |3[ The weak matrix element of the 27 operator seems to vanish in the chiral limit 
with p 2 = 0, while the other elements remain a constant in the limit. These trends in the pion mass 



dependence are reasonably consistent with the prediction of ChPT [|10p. In order to investigate 
the m % and p dependences, we carry out a global fitting for each matrix element for rrr n and p 2 
assuming a simple polynomial form as 



Aoo+Aioml+Aumlp 2 +A m p 2 . 



(3.3) 



The fitting results at zero momentum for each pion mass and at the chiral limit are represented by 
square symbols in the figures. The constant in the fitting result Aoo for the 27 operator is consistent 
with zero within the error, Aqq = —0.028(25), as expected. 

We calculate Re^2 from the weak matrix elements with the Wilson coefficients evaluated by 
NDR scheme, which is converted to RI scheme at the scale p: = 1.44 GeV. The ReA2 obtained 
from the CM and Lab system calculations is shown in Fig. |]. We also plot the results previously 
obtained with H-parity boundary calculation reduction method [p|], and direct calculation with 
ChPT [|12|]. The results strongly depend on the pion mass, and also the momentum. We estimate 
ReA 2 at the physical point, m„ = 0.14 GeV, m K = 0.498 GeV and p = 0.206 GeV, by a global fit 



making the same polynomial assumption as for the case of the weak matrix elements in eq. (3.3). 
The square symbols in the figure denote the fitting results at p = 0.206 GeV. The ReA 2 at the 
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physical point is 2.54(43) x 10~ 8 GeV which is 1.69(28) times larger than the experiment. In order 
to understand the difference from the experiment, we need a more detailed investigation of the 
systematic errors, e.g., finite volume effects, discretization errors, and quenching effect. 
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